Modeling Supply Chains and Business Cycles as Unstable Transport Phenomena 
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Physical concepts developed to describe instabilities in traffic flows can be generalized in a way 
that allows one to understand the well-known instability of supply chains (the so-called "bullwhip 
effect"). That is, small variations in the consumption rate can cause large variations in the production 
rate of companies generating the requested product. Interestingly, the resulting oscillations have 
characteristic frequencies which are considerably lower than the variations in the consumption rate. 
This suggests that instabilities of supply chains may be the reason for the existence of business cycles. 
At the same time, we establish some link to queuing theory and between micro- and macroeconomics. 



Concepts from statistical physics and non-linear dy- 
namics have been very successful in discovering and ex- 
plaining dynamical phenomena in traffic flows [1]. Many 
of these phenomena are based on mechanisms such as 
delayed adaptation to changing conditions and the com- 
petition for limited resources, which are relevant for other 
systems as well. This includes pattern formation such as 
segregation in driven granular media [2] and lane forma- 
tion in colloid physics [3] or biological physics (pedestri- 
ans, ants) [4]. Another example are clogging phenomena 
at bottlenecks in freeway traffic [5] , panicking pedestrian 
crowds [6] , or granular media [7] . In the following study, 
we will focus on the phenomenon of stop-and-go traffic 
[8] and its analogies. 

Recently, economists and traffic scientists have won- 
dered, whether traffic dynamics has also implications for 
the stability and management of supply chains [9,10] or 
for the dynamics of business cycles [12]. To explain busi- 
ness cycles, many theoretical concepts have been sug- 
gested over the decades, such as the Schumpeter clock 
[11]. These are usually based on macroeconomic variables 
such as investment, income, consumption, public expen- 
diture, or the employment rate, and their interactions. In 
contrast, Witt et al. [12] have recently suggested to in- 
terpret business cycles as self-organization phenomenon 
due to a linear instability of production dynamics re- 
lated to stop-and-go waves in traffic or driven many- 
particle systems. In order to illustrate their idea, they 
have transferred a continuous macroscopic traffic model 
and re- interpreted the single terms and variables. 

The author believes that this is a very promising ap- 
proach to understand business cycles, but instead of sim- 
ply transferring macroscopic traffic models, suggests to 
derive equations for business cycles from first principles, 
which means to derive the dynamics on the macroscopic 
level from microscopic interactions. This would also make 
some contribution to the goal of understanding macroe- 
conomics based on microeconomics or, even more, based 
on the "elementary interactions" of individuals. 

In order to make some progress in this direction, we 
will generalize some ideas suggested by Daganzo to de- 
scribe the dynamics of supply chains [9] . Like the work 
by Armbruster et al. [10], his approach is related to traf- 
fic models as well, but he focusses on models in discrete 



space in order to reflect the discreteness of successive 
production steps. In order to be able to reflect the eco- 
nomics of a country, we will have to generalize these ideas 
to complex production and supply networks. This will be 
done in Sec. I. Section II will, then, focus on business cy- 
cles in an sectorally structured economy. This manuscript 
can, of course, only be a first step into the direction of 
simulating complex production processes or the whole 
economy of a country. Some further research directions 
are indicated in the outlook of Sec. III. 

I. MODELLING SUPPLY NETWORKS 

The production units: We will investigate a system with 
u production units (e.g. machines or factories) a,b,c £ 
{1,2,..., u} producing p products i, j, k £ {1,2,... , p}. 
The respective production process is characterized by pa- 
rameters ni and fh In each production step, production 
unit b requires n\ products (educts) j £ {1, . . . ,p} and 
produces fl products i £ {1, . . . ,p}. The number of pro- 
duction steps of production unit b per unit time shall be 
Pb(t). It is the product of the processing (departure) rate 
Pb and the probability pb of the production unit being oc- 
cupied: 

Pb{t) = Hb{t)Pb{pb, c b , s b ) ■ (1) 

Herein, 

p b = Xb/pb (2) 

denotes the channel utilization, A& the feeding rate, Q, 
the number of parallel channels of the production unit, 
and Sb the storage capacity. Changing the number Cb of 
channels or the storage capacity Sb is costly. For the time 
being, we will therefore assume Cb to be constant. The ex- 
pression for Pb{pbi Cb, Sb) is well-known for the stationary 
state of many queuing systems [13]. For example, for a 
M/M/l : (oo/FIFO) process (i.e. one channel with first- 
in-first-out serving, unlimited storage capacity, Poisson- 
distributed arrival times and exponentially distributed 
service intervals), one finds 

p b {pb; 1, oo) = min (p b , 1) , (3) 

so that, in this particular case, we have the relation 
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Pb(t) = fi b (t) mm[p b {t), 1] = min[A 6 (i), fi b {t)] . (4) 

Feeding rates: The feeding rate Xb is not just the prod- 
uct of the concentrations of educts required for a prod- 
uct, as that would be the case for chemical reactions. 
Instead, the feeding rate is determined by the minimum 
arrival rate A 3 ab of required educts j, divided by the num- 
ber n b of educts required for one production step. (The 
difference is comparable to the Probabilistic AND and 
the Fuzzy AND in Fuzzy Logic.) As the arrival rates A J ab 
are given by the input buffer flows v^I^ (where is the 
rate of getting educt j from the the input buffer into the 
production unit b and I b is the number of educts stored 
in the input buffer), the resulting feeding rate is 



Xb(t) = min 

j 



<(t)lj{t) 



(5) 



Note that, due to stochastic variations (which can cause 
queuing effects), efficient production is related to X b < 
r\ib with 1 > r « 0.7 [13]. Therefore, it would be rea- 
sonable to use transport rates v° b according to v b = 
rn b n J b /Il- However, usually there are capacity constraints 
V b in the transport of product j in production unit b, 
which may be reflected by a function U (x, V) such as 



U(x,V) = 



x for x <C V 
V forx>K 



1 + x/V 

Therefore, we will make the specification 



4{t) = u 



ry b (t)n 



(6) 



(7) 



Input and output buffers: Let us assume, each pro- 
duction unit b has input buffers for required educts i 
and output buffers (a warehouse) for the products. We 
will assume the input buffers are filled with Il(t) educts 
i 6 {1, . . . ,p} and the output buffers with B b (t) prod- 
ucts waiting to be delivered. If D l ab (t) denotes the deliv- 
ery flow of products i from production unit a to 6, the 
change of an input buffer stock with time is given by the 
conservation equation 



r ]ji 

-^ = J2 D ib(t)-nip b (t), 



(8) 



as J2 a -^ab(^) 1S t ne number of products i delivered from 
various sources (production units) a and n b Pb is the num- 
ber of educts i used up for production per unit time. 
Analogously, the dynamics of an output buffer stock is 
determined by the equation 



dBi 



dt 



b - = fim)-Y,Di{t) 



(9) 



as f b Pb is the number of newly generated products j, and 
S c -^LW are deliveries to other production units c. The 
delivery flow D l ab is given by the delivery rate v % ab times 
the available products B l a in the buffer: 



DUt) = vib(t)B:(t). 



(10) 



The delivery flows are adapted to the order flows O l ab . 
Ideally, one would have the relation D l ab = O l abl but due 
to capacity constraints V* b , we will again assume 



Kbit) = u 



Ojbit) 
Blit) 



,v: b it) 



(ii) 



which implies D % ab < O l ab . 

Adaptation of capacities: One important aspect of pro- 
duction is the adaptation of production capacities. This 
requires usually a lot of time compared to the other time 
scales of production. The adaptation time 7& of the pro- 
cessing rate /j b is much greater than the other adaptation 
times, so that we will assume the relation 



dub _ J_ 
dt ~ T b 




- Hb 



(12) 



Herein, TL r (y, ...) is some control function reflecting the 
strategy by the production management to adapt the 
processing rate fib, e.g. as a function of the output buffer 
stocks B b or other variables. Normally, it will increase 
with decreasing output buffer stocks B b , but the function 
Wiy, . . .) saturates due to financial, spatial, or technolog- 
ical limitations and inefficiencies in the processing of high 
order flows. In this study, we will assume a function of 
the form 



W(l/z,...) = A 



1 + Bz 



1 + Cz + Dz 2 



(13) 



with z = 1/y and suitable parameters A, B, C, and D. 

Of course, the transportation capacities V b and V* b are 
adapted as well. One may set up separate equations for 
this, but for simplicity, we will here assume the propor- 
tionality relations 

V b j it) = V b ^bit) and V^it) = VMt) . (14) 

Order flows, delivery networks, and price dynamics: 
The flow of orders is basically given by the flow Pbn\ of 
educts required for the production by unit b. If the pro- 
portions q l ab with ^ a q l ab = 1 reflect orders from different 
producers, we have 



o\ b - QlbPbnl . 



(15) 



The fractions q l ab characterize the delivery or supply net- 
work. One can imagine various scenarios. For example, 
the fractions could be modelled by an evolutionary selec- 
tion equation with a selection rate v q [14]: 
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dt 



KM ~ E F a' b (tK> b (t) qU(t) ■ (16) ^ = 1 

dt T b 



It makes sense to relate the fitness F* b to the inverse price 
p l ab of product i, which producer a takes from production 
unit 6: 

Kb® = i/ P i b (t) . (I?) 

There are many pricing strategies, but the common ele- 
ment seems to be that the price changes according to the 
law of supply and demand. Therefore, one conceivable 
specification would be 



Pl b (t)=P 



o QU(t) 



(18) 



where p° is the price when supply D l ab and demand O l ab 
agree. 

Other specifications are, of course, possible as well, as 
the above formulas partly depend on the strategies of the 
human decision makers involved. 



II. MARKETS AND BUSINESS CYCLES 

We will now introduce a simplified and aggregate de- 
scription of an economy. For this, we will summarize a 
whole economic sector by one production unit b and re- 
place input and output buffers by markets for the prod- 
ucts i of the u different economic sectors. To describe 
this, we define the macroscopic stock level of products i 
in the market by 



(19) 



According to the balance equations (8) and (9), we find 
dN t = 

h 



dt 



(20) 



Furthermore, in the arguments of the remaining mathe- 
matical relations, we have to replace input buffer stocks 
l b and output buffer stocks B % b by the respective number 
Ni of products in the market, as the market replaces the 
buffers. In this way, we obtain the feeding rates 



X b (t) = min 



(21) 



which enter P b according to equation (1), the transport 
rates 



rn b (t)n 3 



(22) 



and the processing rates 



p b 



(23) 



For P b « X b , one would have to solve this equation to- 
gether with 



dN, 



^ = £(/j-»j)Mt)niin 



U 



(24) 



For a linear supply chain with n\ = 5 b ,i+i and f b = 5 b ,-, 
(where Sk,i = 1, if k = I, otherwise 0), this leads to 



^=V i (t)N i _ 1 (t)-V i+1 (t)N i (t) 



with 



Vi(t)=m(t)u 



iVi-i(i) 



i-1 



(25) 



(26) 



Interestingly, Eqs. (23) and (25) basically agree with the 
traffic flow model by Hilliges and Weidlich [15], where 
(25) is analogous to the equation for the vehicle den- 
sity and (23) corresponds to the velocity equation. These 
equations are known to behave linearly unstable with re- 
spect to perturbations of Ni(t), if T b exceeds a certain 
threshold which depends on the maximum slope dW/dNi 
[15]. This is the reason why one can frequently observe an 
instability of supply chains, called the "bullwhip effect" 
(e.g. in the "beer distribution game" [16]). A thorough 
analysis shows that the Hilliges- Weidlich model would, 
in fact, be more suitable for the description of supply 
chains than of traffic: First, their velocity equation does 
not contain a convection term, as it should be the case 
for traffic. Second, their equations normally do not show 
metastability [8], i.e. a region in which traffic flow breaks 
down, when a perturbation exceeds a critical amplitude, 
but where smaller perturbations fade away. For this rea- 
son, a supply chain of the above kind is expected to be- 
have either stable or linearly unstable, no matter how 
large the perturbation amplitude is. The most interesting 
point, however, is the reaction of the system to a peri- 
odic perturbation in the consumption rate V p +i(t), as the 
resulting variations in the stock levels are synchronized 
and much slower (see Fig. 1). For this reason, they may 
explain business cycles as self-organized phenomenon. 



III. SUMMARY AND OUTLOOK 

In this contribution, a theory of supply networks has 
been sketched, which may open a new area of econo- 
physics [17]. This theory is developed to help under- 
stand the dynamical phenomena, breakdowns, instabil- 
ities and inefficiencies of production processes and sup- 
ply networks. Future work will have to address questions 
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such as the relevance of the network structure for the re- 
sulting dynamics, possible control strategies, the role of 
the market and pricing mechanism, etc. This specifically 
concerns the choice of the factor r, the specification of the 
functions U, W, and the equations suggested in the para- 
graph on delivery networks and price dynamics. Here, we 
have focussed on an application to markets, i.e. macroe- 
conomics. For this, we have replaced input and output 
buffers by markets for different products [see Eq. (19)]. 
One may view this as some kind of micro-macro link, as 
we have started with equations for single firms (a mi- 
croeconomic level), e.g. Eqs. (8), (9), and ended up with 
equations (20) for market sectors. Assuming a sectoral 
structure of economics, one can relate the resulting equa- 
tions with the Hilliges-Weidlich model, which has origi- 
nally been developed for traffic flow. These equations can 
describe the "bullwhip effect" due to their linear insta- 
bility in a certain regime of operation. The underlying 
mechanism is the slow adaptation of the processing rate 
to changes in the order flows or stock levels in the mar- 
ket. Interestingly, the resulting oscillations in the stock 
levels of the different products i have a characteristic fre- 
quency, which can be much lower than the underlying 
fast variations in the consumption rate. These oscillations 
synchronize among different economic sectors and may 
explain business cycles as a self-organized phenomenon 
with slow dynamics. As the variables in the model are op- 
erational and measurable, the model can be tested and 
calibrated with empirical data. Apart from some equa- 
tions which were not further applied in this study, most 
of the proposed model equations were conservation equa- 
tions, equations given by the product flows, or relations 
derived with stochastic concepts used in queuing theory. 
They reflect the transport and interaction of products, so 
that the physics of driven many-particle systems and of 
complex systems can make some significant contributions 
to the new multidisciplinary field of self-organization phe- 
nomena in production and supply networks. 
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FIG. 1. Variations of the stock levels Ni(t) of u = p = 10 
economic sectors i, triggered by a small perturbation 
5Vu(t) = O.lsin(O.lt) of the constant consumption rate 
V/i = AW (1/20). The model parameters are A = 10, 
B = 0.01, C = 0.01, D = 0.02, 9;- 1 = 0.1, and T b = 90. 
The initial and boundary conditions are Ni(0) = 20 = No(t), 
Vi(0) = V/^W^/N^O)), and Vn(t) = V& + 6Vu(t). 
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